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Abstract. We present the results of a CO survey of 
gravitationally lensed quasars, conducted with the Plateau 
de Bure Interferometer over the last three years. Among 
the 18 objects surveyed, one was detected in CO line emis- 
sion, while six were detected in the continuum at 3mm 
and three in the continuum at 1mm. The low CO detec- 
tion rate may at least in part be due to uncertainties in 
the redshifts derived from quasar broad emission lines. 
The detected CO source, the z = 3.2 radio quiet quasar 
MG0751+2716, is quite strong in the CO(4-3) line and in 
the millimeter/submillimctcr continuum, the latter being 
emission from cool dust. The integrated CO line flux is 
5.96 ± 0.45 Jy-kms -1 , and the total molecular gas mass is 
estimated to be in the range Mh 2 = 1.6 — 3.1 x 10 9 Mq. 

Key words, quasars: general - gravitational lensing - 
quasarsdndividual: MG0751+2716 

1. Introduction 

The measurement of CO line emission in high redshift ob- 
jects has proven to be a fruitful avenue for investigating 
the properties of distant quasars and galaxies. Molecular 
observations address interesting issues such as that of star 
formation in the early universe and how the presence of 
a central massive engine can affect the interstellar ma- 
terial in its host galaxy. Detections of CO have shown 
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that quasars share some properties with luminous infrared 
galaxies, both locally and at high redshift. For example, 
comparison of the molecular, infrared and optical prop- 
erties of the Cloverleaf quasar and the infrared galaxy 
IRAS F10214+4724 (known to harbor a buried quasar 
seen in polarized light) has demonstrated that these two 
objects are nearly identical, except in the optical range 
where the differences can probably be attributed to ob- 
scuration/orientation effects (Barvainis et al. 1995). Such 
findings lend support to theories unifying luminous in- 
frared galaxies and quasars via orientation effects, with 
high redshift infrared galaxies being the luminous counter- 
parts of local Seyfert 2's. However, on physical grounds, it 
also seems likely that IR-selected galaxies and UV-selected 
quasars may differ in their stage of evolution. There are 
currently about 15 well-documented detections of molec- 
ular gas at high redshift (e.g., Combes 2001), of which at 
least 9 are gravitationally lensed systems. 

The advantages of using an intervening "gravitational 
telescope" to boost the fluxes are obvious, with estimated 
magnification factors of up to 100 in the optical. Moreover, 
differential gravitational effects provide an elegant tool 
to probe the size and structure of the molecular mate- 
rial within the quasar. For example, a point-like emitting 
region (rest-frame UV and optical continua from the in- 
ner accretion zone), and an extended dusty molecular re- 
gion (the "torus" ) in the quasar will produce, after grav- 
itational effects from the intervening lens, images with 
different morphologies. Molecular line profiles, reflecting 
intrinsic geometrical and kinematical properties, can be 
particularly useful in understanding the extended struc- 
ture. It should be noted however that a detailed model 
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of the intervening lens must be available to perform the 
transfer from the image plane (observational data) to the 
source plane (intrinsic properties of the quasar). We ap- 
plied this technique for the first time to recover the prop- 
erties of the molecular torus in the Cloverleaf, a quasar 
at z = 2.56 (Kncib et al. 1998), comparing HST images 
and IRAM interferometer CO maps (Alloin et al. 1997). 
The CO-emitting region in the quasar was found to be a 
disk or ring-like structure orbiting the central engine at 
a radius between 75 and 100 pc, with Keplerian veloc- 
ity around 100 km s _1 . The effective resolution resulting 
from this technique turned out in this case to be about 20 
times smaller than the synthesized beam size of the CO 
interferometer data. 

Such significant benefits - flux boosting and increased 
effective angular resolution - have led us to focus our at- 
tention on gravitationally lensed systems and to conduct 
a CO survey of these objects. Another possible benefit of 
lensed versus unlensed objects is the potential for higher 
flux boosting for higher-J transitions (differential mag- 
nification). Since experience has shown that the best se- 
lection criterion for CO detection is the presence of de- 
tectable far-IR or submm/mm dust continuum emission, 
we started with a continuum survey of the known lensed 
quasars using the IRAM 30m radio telescope and the 
JCMT (Barvainis & Ivison 2002). A high dust continuum 
detection rate encouraged us to pursue a CO search with 
the IRAM interferometer (Guilloteau et al. 1992). 

Since this project was started early in the continuum 
survey, we were not at that time able to make a general 
selection based upon submillimcter flux. Instead, the sam- 
ple consisted of most of the then-known lensed quasars 
having optical redshifts measured to good accuracy. We 
observed 18 gravitationally lensed quasars, with redshifts 
in the range 1.375 — 3.595. 

However, reliable systemic redshifts remain a major 
difficulty for CO searches at high z because currently 
available redshifts are mostly derived from highly ion- 
ized species in the quasar broad line region. As this re- 
gion is often coupled to a high velocity wind, redshifts 
derived this way have been found to be blucshiftcd up 
to 1200 km s _1 with respect to the systemic velocity of 
the host galaxy and the molecular environment of the 
quasar probed by CO measurements. A typical offset is 
600 km s _1 , but there is wide dispersion from one ob- 
ject to another. Meanwhile, spectrometer bandwidths in 
the millimeter domain are too narrow (~ 1500 km s^ 1 at 
3mm) to fully span this redshift uncertainty using a single 
central frequency setting. The combination of these two 
facts makes it likely that some CO lines will be missed 
in the course of a survey. In the case of the present sur- 
vey, whenever the quasar redshift was from highly ionized 



species we applied a 600 km s _1 redshift increment to 
search for its CO emission. We are fully aware that this 
offset, although statistically meaningful, may be just in- 
correct for some individual quasars. 

In Sect. 2 we describe the sample of gravitationally 
lensed quasars and the acquisition and reduction of the 
interferometer data set. Results, both in CO line emission 
and in the 1mm and 3mm continua, are also presented in 
Sect. 2 for the entire sample. In Sect. 3, we discuss the 
general results of the CO survey, and in Sect. 4 consider 
the detection of MG0751+2716 in the CO(4-3) transition 
in more detail. Conclusions and future prospects are given 
in Sect. 5. 



2. The sample: acquisition and reduction of the 
CO interferometric data set 

The gravitationally lensed quasar sample is presented in 
Table 1. The coordinates generally refer to the brightest 
quasar image in the optical. Redshifts have been corrected 
in some cases using the technique described in the previous 
section, except for MG0751+2716 where an initial detec- 
tion allowed refinement of the value to the center of the 
line in followup observations. The centering frequencies 
given for the 3mm and 1mm windows correspond to var- 
ious CO transitions from CO(2— 1) to CO(9— 8), depend- 
ing on the quasar redshift and the window considered: the 
targetted CO transitions in the 3mm window are specified 
for each source in Table 1. The "seeing" estimates corre- 
spond to the 3mm data set. Some targets were observed 
on several different observing runs; seeing values, numbers 
of antennas, and hours spent per exposure are provided. 
The observed spectral bandwidth was either 560 MHz or 
595 MHz. 

Table 2 lists the CO transition observed (3mm), the 
"channel rms" (rms/beam/100 km s^ 1 , in mJy), the line 
flux (for MG0751+2716), and the 3mm and 1mm contin- 
uum results. Sufficient sensitivity and bandwidth for line 
measurements were only available at 3mm. Table 2 also 
lists 850 /im continuum fluxes or upper limits derived from 
SCUBA observations at the JCMT (Barvainis & Ivison 
2002). Note that B1600+434 and B1030+074 are strong 
radio loud quasars whose millimeter/submillimeter con- 
tinua are consistent with pure synchrotron radiation. The 
other submillimeter-detected sources are likely to be dom- 
inated by dust emission. The continuum flux ratio between 
3mm and 1mm for MG0751+2716 is not typical of dust 
emission, though the ratio between 1mm and 850 /<m is. 
This is probably caused by some residual steep-spectrum 
synchrotron emission contributing to the 3mm (and to a 
lesser extent the 1mm) flux. 
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Object 


RA(2000) 


Dec(2000) 


z 


I^3 mm 


Seeing 


Beam sizc/PA 




Beam sizc/PA 


Antennas 


Exposure 










(GHz) 


(") 


(" x "/ a ) 


(GHz) 


(" x "/") 




(Hrs) 


0047-2808 


00:49:41.87 


— 27:52:25.7 


3.595 


100.335 


1.6:1.0 


14.1 x 3.5/11 


225.661 




4;5 


3.5;1.6 


UM673 


01:45:17.22 


— 09-45. 1 9 3 


2.730 


92.706 


0.75 


A r i v 1 Q / n 


216.260 


"1 Q v 7 / 1 QQ 


5 


4.2 


MG0751+2716 


07:51:41.46 


+27:16:31.4 


3.200 


109.778 


1.0;1.1;0.6 


9.4 x 5.6/32 


246.898 


3.1 X 2.1/14 


4;3;4 


3.4;3.6;2.0 


SBS0909+523 


09:13:00.76 


+52:59:31.5 


1.375 


97.068 


0.2 


2.4 x 1.7/155 


242.639 


0.9 X 0.7/147 


5 


4.7 


RXJ0911+055 


09:11:27.50 


+05:50:52.0 


2.807 


90.831 


0.7 


6.5 X 6.0/53 


211.886 


3.0 X 2.0/127 


4 


6.3 


Q1009-0252 


10:12:16.09 


-03:07:03.0 


2.746 


92.311 


0.7 


13.0 X 6.2/9 


215.337 


5.5 x 2.9/15 


4 


3.1 


J0313 


10:17:24.13 


-20:47:00.4 


2.552 


97.352 


0.9 


17.4 x 6.3/7 


227.097 


8.0 x 2.4/7 


4 


4.0 


B1030+074 


10:33:34.0 


+ 07:11:26.1 


1.535 


90.942 


1.1;2.3 


9.9 X 5.8/49 


227.325 




5;5 


1.1;3.2 


HE1104-1805 


11:06:33.45 


-18:21:24.2 


2.326 


103.967 


1.7 


13.5 X 4.5/8 


242.529 




3 


3.2 


PG1115+080 


11:18:16.96 


+07:45:59.3 


1.723 


84.663 


0.7 


7.8 X 5.5/53 


211.630 


2.8 x 2.2/60 


4 


5.8 


1208+1011 


12:10:57.16 


+09:54:25.6 


3.831 


95.434 


0.5 


2.5 X 1.4/23 


214.637 


1.1 x 0.6/24 


4 


6.3 


HST14176+530 


14:17:36.51 


+52:26:40.4 


3.403 


104.711 


0.8|2.0 


6.6 X 4.8/50 


235.501 


2.4 x 1.7/98 


4;5 


3.7;5.9 


SBS1520+530 


15:21:44.83 


+52:54:48.6 


1.860 


80.608 


1.1;0.7;0.7;1.2 


10 X 5.7/123 


241.773 


5.3 X 2.5/63 


4;3;5;5 


4.4;3.8;0.33;2.3 


B1600+434 


16:01:40.45 


+43:16:47.8 


1.589 


89.045 


0.2;1.4;2.0 


8.3 X 7.6/138 


222.583 


1.0 X 0.6/46 


5;4;4 


1.9;6.1;2.9 


B1608+656 


16:09:13.96 


+65:32:29.0 


1.394 


96.298 


0.4 


3.0 x 1.7/83 


240.713 


1.2 X 0.7/75 


5 


3.9 


2016+112 


20:19:18.15 


+ 11:27:08.3 


3.282 


107.670 


0.7;0.6|1.3 


7.8 X 5.0/45 


242.156 


3.0 X 2.4/21 


5;5;4 


4.9;2.6;4.9 


HE2149-2745 


21:52:07.44 


-27:31:50.2 


2.033 


114.011 


3.6 


24.1 x 4.6/7 


227.938 




5 


3.0 


Q2237+0305 


22:40:30.14 


+ 03:21:31.0 


1.696 


85.511 


1.0 


6.5 X 4.1/2 


213.749 


2.6 X 1.6/176 


4 


5.1 



Table 2. Results of the CO survey 



Object 


Target line 


Line rms 


Line flux 


3mm Scont 


1mm Scont 


850^m Scont 




(mjy/beam) 


(mJy/lOOkms- 1 ) 


(Jy-kms -1 ) 


(mjy/beam) 


(mjy/beam) 


(mJy) 


0047-2808 


CO(4-3) 


3.3 




0.9 ± 0.9 


na 


< 7.0 


UM673 


CO(3-2) 


2.0 




0.7 ± 0.5 


na 


12.0 ±2.2 


MG0751+2716 


CO(4-3) 


1.9 


5.96 + 0.45 


4.1 ± 0.5 


6.7 ± 1.3 


25.8 ± 1.3 


SBS0909+523 


CO(2-l) 


1.1 




0.5 ± 0.3 


0.0 ± 1.3 


< 5.5 


RXJ0911+055 


CO(3-2) 


1.2 




1.7 ± 0.3 


10.2 ± 1.8 


26.7 ± 1.4 


Q1009-0252 


CO(3-2) 


1.8 




0.1 ± 0.5 


0.3 ± 4.0 


na 


J0313 


CO(3-2) 


1.8 




0.0 ± 0.5 


-6 ± 4 


na 


B1030+074 


CO(2-l) 


3.3 




184 ± 1.9 


na 


na 


HE1104-1805 


CO(3-2) 


6.8 




1.9 ± 2.0 


na 


14.8 ±3.4 


PG1115+080 


CO(2-l) 


1.3 




-0.3 ± 0.3 


-0.9 ± 1.2 


3.7 ± 1.3 


1208+1011 


CO(4-3) 


1.3 




0.2 ± 0.3 


4.2 ± 1.9 


8.1 ±2.0 


HST14176+530 


CO(4-3) 


1.2 




0.1 ± 0.3 


-5 ± 3 


< 3.5 


SBS1520+530 


CO(2-l) 


1.5 


hint 


0.0 ± 0.2 


-0.5 ± 1.4 


9.4 ±2.6 


B1600+434 


CO(2-l) 


1.1 


hint 


25 ± 0.3 


12.6 ± 2.3 


7.3 ± 1.8 


B1608+656 


CO(2-l) 


1.5 


hint 


8.1 ± 0.4 


0.8 ± 2.2 


8.1 ± 1.7 


2016+112 


CO(4-3) 


0.8 




1.8 ± 0.2 


1.1 ± 1.0 


< 4.8 


HE2149-2745 


CO(3-2) 


10. 




-2.8 ± 2.7 


na 


8.0 ± 1.9 


Q2237+0305 


CO(2-l) 


1.4 




0.1 ± 0.35 


-1.9 ± 2.1 


3.9 ± 1.2 



All of the continuum detections, and the one line de- 
tection, are within <~ 1" of the optical positions listed 
in Table 1, except for two cases. For 2016±112, the off- 
set is <5RA = — 1'/5 and SDec = 3" 5 from the observed 
coordinates. However, checking the NASA Extragalactic 
Database (NED) we found improved coordinates which 
are within about 1" of the continuum source in both RA 
and Dec. In the case of RXJ0911±055, the offset of the 
continuum source (a mean of the 3mm and 1mm contin- 
uum sources) is <5RA = 2'.'0 and SDec = 2'/l. The origin 
of this offset, which is larger than typical optical position 
errors, is unknown at present. 



3. Results of the CO survey 

Out of 18 sources observed, we obtained: (a) one strong 
CO line emitter, MG0751±2716, detected in CO(4-3), (b) 
three marginal detections (to be investigated further) and 
(c) 14 non-detections. Generally speaking, strong submm 
emission from dust is a good predictor of strong CO emis- 
sion. Of the two strongest dust sources, RXJ0911±055 
(S850 = 26.7 mJy) and MG0751±2716 (S 85 o = 25.8 mJy), 
only the latter was detected in CO. The line, with a peak 
flux of 18 mJy, is one of the strongest known among high-z 
sources. Seven sources had either no measurements or only 
upper limits at 850 /xm, but 6 additional sources had mod- 
erately strong submm continuum detections in the 8 — 15 
mJy range. 
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We believe some CO lines, such as the strong one 
expected from RXJ0911+055, and some among the 
intermediate-strength submm sources, may have been 
missed because of the uncertainty in the redshift of the 
molecular lines combined with the narrow observing band- 
width. In fact, during our 1999 first search for CO in 
MG0751+2716, the CO(4-3) transition was detected on 
the edge of the bandpass, at a velocity offset of —600 
km s _1 with respect to our original guess for molecular 
emission. Thanks to the strength of this line we were able 
to identify its presence and obtain new observations at 
the appropriate frequency, which fully confirmed the 1999 
measurement. However, this would not work in the case of 
fainter CO line emitters. This difficulty will only be fully 
resolved by using broadband backends in the future. In 
the meantime, we plan to rcobserve at least RXJ0911+055 
with flanking bandpasses to cover more redshift space. 

Continuum emission was detected at 3mm in 6 targets, 
with 3 of those also detected at 1mm (one being particu- 
larly strong, the radio loud quasar B1030+074). 

4. The CO line in MG0751+2716 

A strong line in the CO(4— 3) transition was detected in 
MG0751+2761. A set of the channel maps, with velocity 
steps of 100 km s _1 , is shown in Fig. 1 for the CO(4— 3) 
transition. It convincingly reveals the CO(4— 3) emitting 
region at a location very close to the quasar optical coor- 
dinates (offset by SRA = 076, <5Dec = 075). The CO(4-3) 
line profile is displayed in Fig. 2. A Gaussian fit provides 
the following parameters: FWHM of 390 ±38 km s _1 , peak 
frequency at 109.778 ±0.005 GHz (corresponding to a red- 
shift of 3.200), and intensity of 5.96 ±0.45 Jy-km s" 1 . The 
CO (9— 8) transition was not detected in the 1mm window. 
A Gaussian fit with fixed line width and central position 
analogous to the CO(4— 3) transition showed only contin- 
uum at the position. A local peak close to 100 km s _1 in 
the object's velocity frame is below the 3cr level. By com- 
bining both sidebands of the 1mm receiver, we obtain a 
continuum detection of 6.75 ± 1.32 mJy. 

In order to compute the CO line luminosity and to- 
tal gas mass, a lensing correction must be applied. In 
spite of several efforts to model the lens system towards 
MG0751+2716 (Lehar et al. 1997; Tonry & Kochanck 
1999), additional work remains to be done: the lens ap- 
pears to be a quite complex system which requires more 
shear than accounted for by the lensing galaxy (identi- 
fied as G3 in Lehar et al. 1997). This extra shear might 
come from what appears to be a group or cluster of galax- 
ies, indicated by the large number of galaxies in the field 
around the quasar. Very recent modeling by J. Lehar 
and B. McLeod (2001, private communication), based 
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Fig. 1. Channel maps of the CO (4 — 3) line emission in 
MG0751±2716 with a contour spacing of 2.5 mJy/beam 
(the zero contour is not shown). The cross indicates the 
optical position of the quasar. Two different spectral se- 
tups were merged which overlap in the [—75,725] km s _1 
channels, resulting in three different rms levels over the 
spectral range (see Fig. 2 error bars). For the deconvolu- 
tion of each channel, the appropriate combined UV cover- 
age was taken into account. The continuum flux has been 
subtracted. 

on HST optical imaging, provides an estimated optical 
magnification of 16.6. While differences between optical 
and CO /submillimeter net magnifications can be expected 
in some cases because of the difference in source sizes 
(subpc versus tens to hundreds of pc), models suggest 
that for optical magnifications less than about 20 the dif- 
ferences are generally not expected to be large (see Fig. 
1 of Barvainis and Ivison 2002, and associated discus- 
sion). Therefore, we correct the CO emission for a mag- 
nification of factor of 16.6 and derive a line luminosity 
L' CQ = 3.9 x 10 9 K-kms- 1 pc 2 (H = 75 km s" 1 Mpc" 1 , 
?o-0.5). 

The molecular material emitting in the CO (4— 3) tran- 
sition is most likely close to the quasar, as in the case of 
the CloverLeaf (Kneib et al. 1998) and APM 08279+5255 
(Downes et al. 1999), and it is also probably dense and 
warm. In order to calculate the total molecular gas mass 
we have considered two values for the ratio M H2 /£' co : 
0.8Mq (K-kms -1 pc 2 ) -1 , following Downes & Solomon 
(1998) for nuclear rings in ultraluminous galaxies, and 
O.4M (K-kms- 1 pc 2 )" 1 following Barvainis et al. (1997) 
for the molecular "torus" in the Clovcrlcaf. These conver- 
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Fig. 2. Spectra in CO(4 - 3) (upper panel) and CO(9 - 8) 
(lower panel) for MG0751+2716, binned into 100 km s" 1 
wide channels. 

sion factors yield a total molecular gas mass in the range 
M H2 = 1.6-3.1 X 1O 9 M . 

It is of interest as well to compare the gas mass to 
the dynamical mass derived from the observed CO line 
width. With a derived Keplerian velocity of 400/sin(i) km 
s _1 , and assuming that the molecular gas is located at a 
radius of 200 pc (see above references for the CloverLeaf 
and APM 08279+5255), we obtain a lower limit of M dyn > 
1.8 x 10 9 Mq, consistent with the derived value of Mh 2 . 



Interferometer, all of the observations related to this project. 
We also thank Ski Antonucci for important early contribu- 
tions. Joseph Lehar and Brian McLeod generously provided 
an estimate of the magnification of MG0751+2716 in advance 
of publication. DA and RB wish to thank the IRAM insti- 
tute in Grenoble for hospitality during visits there. The NASA 
Extragala ctic Database (NED) and the CASTL eS compilation 
of lenses ( jittp:/ /cfa- www. harvard. edu/castles[ ) were used ex- 
tensively in the course of this work. 
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5. Concluding remarks and future prospects 

Though the present survey yielded a low detection rate 
in CO, there are several new lensed quasar candidates 
yet to be observed based on their strong 850 /urn con- 
tinua, recently discovered in the course of the submillime- 
ter survey by Barvainis & Ivison (2002). Supplementary, 
expanded-frequency observations of some sources (most 
notably RXJ0911+055) may turn up more CO detections 
from the present source list. 

As for MG0751+2716, the centimeter radio source has 
four components connected by arcs (Lehar et al. 1997), 
and in the optical it appears as a 1" diameter Einstein 
Ring. A primary driver for this project was to find lensed 
sources that could be spatially resolved in CO line emis- 
sion. This is currently possible for MG0751+2716 using 
the PdBI, and, like the Cloverleaf, reconstruction of the 
molecular source structure and kinematics on very small 
angular scales using the lensing properties may prove to 
be quite interesting. 
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